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0^ Abstract. Wc calculate the renormalized vacuum polarization and stress tensor for 

''"^ a massless, arbitrarily coupled scalar field in the Hartle-Hawking vacuum state on the 

^ horizon of a Schwarzschild black hole threaded by an infinte straight cosmic string. This 

^ ^ calculation relies on a generalized Heine identity for non-intcgcr Lcgcndrc functions 

I which we derive without using specific properties of the Legendre functions themselves. 
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1. Introduction 

There is a long and successful history of calculating the renormalized vacuum 
polarization and stress-energy tensor on spherically-symmetric black holes (see P, El 
El m El El H IB]), however, in the astrophysically significant case of the Kerr- Newman 
black hole, such calculations have proved intractable. This is due to the fact that 
all the calculations in the spherically symmetric case have relied heavily on specific 
properties of the Legendre functions (mainly the Legendre Addition Theorem and the 
Heine Identity) and also the applicability of the Watson-Sommerfeld Formula to obtain 
a mode-sum expression that is numerically tractable. In the axially symmetric case, 
however, addition theorems for the angular functions do not exist nor do we have an 
analog of the Heine identity for the angular functions that arise. We will show, in the 
framework of the Schwarzschild black hole threaded by a cosmic string, that we can 
proceed in the axially symmetric case by obtaining useful summation formulae based on 
the Hadamard structure of the Green function without using specific properties of the 
mode-functions themselves. In this paper we will adopt this approach in order to derive 
a generahzation of the standard Heine identity, adapted to the non-integer Legendre 
functions which arise in the separation of the wave equation on the background metric 
of a Schwarzschild black hole threaded by a cosmic string. Furthermore, we apply this 
formula to the calculation of the vacuum polarization and renormalized stress tensor for 
a massless scalar field on the horizon of the black hole. 

2. The Green's Function 

The cosmic string is modeled by introducing an azimuthal deficit parameter, a, into the 
standard Schwarzschild metric. We may describe this in coordinates (t, r, 9, 0), where 
is periodic with period so we may take G [0,27ra), in which the line element is 
given by 

ds^ = -(1 - 2M/r)dt^ + (1 - 2M/r)"Mr2 + rM^^ + sin^ (1) 

For GUT scale cosmic strings a = 1 — 4yU where ^ is the mass per unit length of the 
string so we shall assume < a < 1. We can alternatively define a new azimuthal 
coordinate by 

= a0 (2) 

so that is periodic with period 27r and we may take G [0,27r). In coordinates 
(t, r, 6', 0), the line element is given by 

ds^ = -(1 - 2M/r)dt^ + (1 - 2M/r)-Mr2 + rM^^ ^ ^2^2 ^^^^2 q^^2 _ 

We shall consider a massless, minimally coupled scalar field, (p, in the Hartle- 
Hawking vacuum state. Since this is a thermal state, it is convenient to work with 
the Euclidean Green's function, performing a Wick rotation of the temporal coordinate 
t ^ —IT and eliminating the conical singularity at r = 2M by making r periodic with 
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period 27r/K where /t = 1/(4M) is the surface gravity of the black hole. The massless, 
minimally coupled scalar field, satisfies the homogenous wave-equation 

□¥.(r,r,^,0) = O, (4) 
which can be solved by a separation of variables by writing 

(^(r, r, e, (j)) ~ e'"""+'"^'^P(^)i?(r) (5) 
where P{0) is regular and satisfies 

while R{r) satisfies 

_ 2Mr)^ - A(A + 1) - I^(^) = 0. (7) 

The A(A + 1) term arises as the separation constant. The choice of A is arbitrary for 
to satisfy the wave equation but requires a specific choice in order for the mode- 
function P{9) to satisfy the boundary conditions of regularity on the poles. In the 
Schwarzschild case (no cosmic string, a = 1), regularity on the poles means that X = I, 
i.e the separation constant is /(/ -|- 1). In the cosmic string case, the appropriate choice 
of A that guarantees regularity of the angular functions on the poles is 

A = / — |m| + \m\/a. (8) 

(We note here the dependence of A on / and m. However, to avoid being typographically 
cumbersome, we will leave it implicit in the remainder of this paper.) With this choice 
of A, the angular function is the Legendre function of both non-integer order and non- 
integer degree, 

P(^) = p;H/«(^Qg^)^ 

We have shown in the Appendix that these angular functions satisfy the following 
normalization condition, 

^r-K^i/ Jcos ^)P^7:/:|^|/ Jcos e)d{cos 6) = 

2 r(A-|m|/a + l) 
(2A + l)r(A+|m|/a + !)"'■ ^ ' 

The periodicity of the Green's function with respect to (r — r') and {(f) — (j)') with 
periodicity 2n/K and 2n, respectively, combined with Eq.(lO) allow us to write the 
mode-sum expression for the Green's function as 

^^^'^^-8vr2 2^ ^ 2^ 2^^^^ + ^^r(A-|m|/« + l) 

n=— oo m=— oo l=\m\ 

p-l'"l/"(cose)P,-l'"l/"(cose')XnA(r,r'), (11) 
where Xnx{f^,f^') satisfies the inhomogeneous equation. 



. dr 



[r' _ 2Mr)- - A(A + 1) - , \xnx{r, r') = —5{r - r'). (12) 

dr — 2Mr ) a 
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It is convenient to write the radial equation in terms of a new radial variable 
?7 = r/M — 1, the radial equation then reads 

- '^^) - ^ 1) - i^i^-C'"') - -^^(" - "')■ (^«) 

where we have used the fact that k = 1/4M. For n = 0, the two solutions of the 
homogeneous equation are the Legendre functions of the first and second kind. For 
n 7^ 0, the homogeneous equation cannot be solved in terms of known functions and 
must be solved numerically. We denote the two solutions that are regular on the horizon 
and infinity (or some outer boundary) hy p^^\r)) and q^xKv)^ respectively. A near-horizon 
Frobenius analysis for n ^ shows that the indicial exponent is ±|?t,|/2, and so we have 
the following asymptotic forms: 

pi'Kv) - iv - ip^' v^i, 

Defining the normalizations by these asymptotic forms and using the Wronskian 
conditions one can obtain the appropriate normalization of the Green's function: 

^Px{v<)Qx{v>) n = 0, 

Xnxiv,v') = ^ " 1 I , I , (15) 

-Px{v<)Qx{v>) n=^0. 



2\n\aM 



3. Renormalization 



From the asymptotic forms (14), it is clear that taking one point on the horizon (with 
the other outside) means that all the n ^ modes vanish. The unrenormalized Green's 
function then reduces to 

Gix x') = ^ V e™(^-*') V (2A + n HA + H/g + l) 

^^"^'"^^ 327r2M2« 2^ ^ Z^^^^ + ^^r(A-|m|/« + l) 

m=-oo l=\m\ VII/ / 

P;:l™l/"(cos^^)P^"l™l/°(cos^^')QA(r/). (16) 

In the absence of the cosmic string, i.e. a = 1, this sum can be done using a combination 
of the Legendre Addition Theorem and the Heine Identity, yielding a closed form 
expression [1], that, of course, diverges in the coincidence limit. 

To renormalize we follow the spirit of the Christensen-DeWitt [9] point-separation 
method subtracting the geometrical singular part of the Green's function for small 
separations in the coordinates. This method relies on the universal Hadamard 
singularity structure of the Green's function which ensures that once we have subtracted 
the geometric singular part, 'f//cr + \^lncr', the coincidence limit of the remainder W 
is finite. To calculate {(f'^)ren in the massless Ricci-flat case, the only subtraction term 
needed is 
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where s is the geodesic distance between x and x'. Having aheady taken one point {x' 
say) on the horizon we choose to separate radially, placing x at 2M + e; the geodesic 
distance is then given by 



2MH 



dr' 



,3 5,4 



,1/2 



2+ ' ' 



2 



^448M3 9216M4 ^ 45056M5 ^ ^ 



6M 80M2 



Note that the leading term behaves like (Ax)^/^ rather that the usual Ax since the 
metric is singular at r = 2M so g^^ = 0(e), but this is simply a coordinate singularity 
and in no way affects the validity of the geometrical subtraction. The subtraction terms 
up to 0(e) are 

^^^"^ " 327r2Me ~ 19M^ + " 327r2M2(r/ - 1) ~ 19MW^ + ^^"^ ~ ^^^^ 



Taking the partial coincidence limit (p ^ (f)' , 9 ^ 9' in (16), subtracting the 



renormalization terms (19) and taking the limit as — t- 1 (e — 0) gives us the 
renormalized expression for the vacuum polarization. However, we now face our 
fundamental challenge as the limit cannot be taken in a meaningful way in the forms 
given. In order to take the limit, we must either invent a way to write the divergent 
terms as an appropriate mode-sum and do a mode-by-mode subtraction or alternatively. 



we must attempt to perform the sum in (34) so that we can write the Green's function 



in closed form. In [TT] we follow the former route to obtain numerical values of (<^^)ren 
outside our black hole. However, when possible, the latter is clearly preferred since it 
can give us a simple closed-form solution for the renormalized vacuum polarization. In 



the next section, we derive a formula that will allow us to perform the sum (16). 



4. Generalized Heine Identity 

We recall that for real arguments, the Heine Identity for the Legendre Functions is 

oo ^ 

^(2/ + l)Pz(*)gKC) = J^^y (20) 

valid for |\E'| < \(\. Writing = C0S7 = cos^cos^' + sin sin cos (0 — 0'), and using 
the Legendre Addition Theorem, 

PKCOS7) = e''^^'^-'''^jf^,Pr{^os9)Pr{cos9% (21) 

we can re-write the Heine Identity as 

E E e-^^-^\2l + l)|^Pr(cos^)Pr(cos^')QKC) = jy^^y (22) 
1=0 m=-i V "T vs ij 

Below we shall generalize this form of the Heine Identity to the non-integer Legendre 
functions that arise in the mode-sum of a cosmic string space-time. It is this form 
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rather than a generahzation of (20) that is useful to us since the fundamentaly axially- 
symmetric nature of the Green's function for a ^ 1 means that its radial part will 
always depend on m (through A), and therefore the m-sum over the angular functions 
can never be done independently of the radial part. 

In the Appendix, we have derived other identities involving Legendre Functions of 
non-integer order and/or degree. We have done so since, although they are not directly 
useful here, they turn out to be useful in other contexts, for example, determining the 
Hadamard singularity structure of a 3D Green's function on a cosmic string black hole 
space-time [TT] . 



Returning to the problem at hand, we want to construct a generalization of (22). 
In addition, we would prefer to do this in such a way that the derivation does not rely 
on specific properties of the Legendre functions themselves, since ultimately we wish 
to apply the same techniques in Kerr space-time. There is a very natural way to do 
this by equating equivalent expressions for the same Green's function on a conveniently 
chosen space-time. In our case, the most appropriate case to consider is that of an 
infinite cosmic string in otherwise flat space-time. Of course, normally in this case one 
would exploit translational invariance along the string and work in cylindrical polar 
coordinates but for our purposes we choose to work in spherical polar coordinates, so 

ds' = dr^ + dr^ + r^de^ + aV^ sin^ Odcj)^, (23) 

where is periodic with period 2tt. The Green's function for a massless, minimally 
coupled scalar field on this space-time has a closed-form solution, which can be obtained, 
for example, by separating in cylindrical polar coordinates 



where 



G(x, x') = — ^^^Hx/a) ^ .24) 

' Stt^q; pp' sinh x(cosh(x/ a) — cos(0 — 0')) ' 

(r — r')^ + r'^ — 2rr'cos6'cos6'' 

coshx = 25 

2pp' 

and p = rsin^ is the cylindrical polar coordinate. 

One can also obtain the spherical polar mode-sum expression for this Green's 
function: 

4vr2a7o ^f^^ r(A - |m|/a + 1) 

p-l"l/'^(cos0)p-l"^l/"(cos^^')^^A(r,r'), (26) 
where Ri_jx{r,r') is the ID Green's function satisfying 

i i^'i) - - + = - ^^^^ 

The solution of this equation can be given in terms of Modified Bessel Functions, 

^-A(r,r) = , (28) 
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where / and K are the Modified Bessel Functions of the first and second kinds, 
respectively. The integral over u can now be performed using the identity [16] 



cosw(r - T')Ix+i/2iujr<)Kx+i/2{(^ry)duj = \i/2 <?a(C) (29) 
valid for A > — 1, where Qx is the Legendre function of the second kind and 

C = L_ . (30) 

Thus, the spherical polar mode-sum expression for the Green's Function may be written 
as 

m=-oo l=\m\ \ I II ; 

P;rl'"l/"(C0S^)P;:I'"I/"(C0S^^')QA(C)- (31) 

We now equate this Green's Function with the equivalent closed form expression 
(|24|) to obtain the following generalized Heine Identity: 



^m,,-,) ^ (2A + 1) [;;^ + + || p,-'"^'"(cOsg)P,-|-l/"(cOsgOgA(C) 

sinh(x/a;) ^^^^ 



sin^sin^'sinhx(cosh(x/a) — cos(0 — 0')) 
where 

, C ^ cos 6* cos 6*' ^ , , ,, ■ n ■ ni /oo\ 

coshx= — = COS7 + (coshx — cos(0 — jj smt^ smc^ . [66) 

sin d sin t'' 



This identity, valid for > | cos7|, is completely analagous to (22) and reduces to it in 
the a — )• 1 limit. 

5. Calculation of {(p^)ren 

Returning now to our calculation of the vacuum polarization on the horizon of the 
Schwarzschild black hole threaded by a cosmic string, we see that our generalized Heine 



Identity (32) has precisely the right form to enable us to perform the mode-sum in the 



Green's function expression ( 11 ). The Green's function, where one point has been taken 



to lie on the horizon, may now be written as 

Gix x') = ^ ^^Hx/a) .34. 

^' ^ 327r2M2asin^sin0'sinhx(cosh(x/a) -cos(0-0')) 
For purely radial separation, we take (j)' cj), 9' 9, to yield 

1 sinh(x/a 
327r2M2a; sin^ 6'sinhx(cosh(x/a) — 1 
where we now have 

coshx = l + ^^^^. (36) 
sm 9 
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From Eqs.(|19|) and (|35j), the renormahzed vacuum polarization is given by 

G{r,e;2M,e)-Gs^n9{r;2M) 



2 \ horizon 
ren 



hm 



hm 



sinh(x/a) 



1 



-327r2M2asin^6'sinhx(cosh(x/a) - 1) 
It is now convenient first expand in terms of x 



327r2M2(r^ - 1) 1927r2M2 



sinh(x/a) 



then 



so we obtain 



sinh x(cosh(x/a;) — 1) 
(r/-l) 



2a 



+ 



1 

6a 



cosh X = 1 + 



sin^^ 



sin^^ 



2{r] 



sinh(x/a) 



sinh x(cosh(x/a) — 1) 



a sin^ ^ _j_ -'- I 
{rj — 1) 6a 



1) 12 ^ ' ' 



1 - a^) + 0(r/ - 1)^ 



(37) 



(38) 



(39) 



(40) 



Substituting (40) into (37) and taking the limit, we arrive at the renormalized vacuum 



polarization on the horizon, 



■x2\horizon 
'ren 



1 



1 



a 



192ti^M^V ' a^sinHJ' 
The first term here is simply the Candelas result [1], which, of course we recover 
when a — )■ 1, and the second term represents the contribution due to the presence of the 



cosmic string. The result (41) has in fact been calculated by Davies and Sahni [10] but 
for a restricted set of a values such that (1/a) is an integer. For this restricted class of 
values, the cosmic string solution is obtained as a sum of (1/a) images of the solution in 
the absence of a string. The method of images cannot be applied for general azimuthal 
deficits, however, so our method presents a generalization and independent derivation 
of the result of Davies and Sahni. 

We have plotted the vacuum polarization in Figure [T] for a range of values of a. We 
see that the presence of the string increases {'p'^)^°n^°^ everywhere. In the equatorial 
plane it is increased by a factor of l/a^ from its value in the absence of a string. As we 
approach the axis (<,3^)j!e"^°"' diverges in a non-integrable manner as 

1 (l-a^) 1 



2 \ horizon 



487r2 



a^ 



(2M sin( 



sin^ 0, 



(42) 



as one would expect a distance (2Msin^) from a flat-space cosmic string. The region 
over which the string dominates may be characterised by the range of cos 6 (near ^ = 0) 



at which (<^^)re"^°"' exceeds twice its value in the equatorial plane, this is given by 



1 — (cos I 



(2-a2)i/2 



(l-a) + 0((l-a)2)). 



(43) 
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Figure 1. {<p'^)ren^°" (ill units M = 1) as a function of cos 6 for a range of values of a. 
For a ^ I, {'p'^)ren^°'^ diverges at the poles since there is a curvature singularity there 
due to the presence of the string. increases everywhere as a decreases. 



6. The Stress Tensor for a Massless Scalar Field with Arbitrary Coupling 



In this section, we shall employ our generalized Heine identity (20) in order to calculate 
the renormalized expectation value of the stress energy tensor for a massless scalar field 
with arbitrarily coupling to the Ricci scalar in the Hartle-Hawking vacuum state. This 
quantity is of particular physical interest in obtaining the back-reaction on the spacetime 
geometry in the semi-classical approximation to quantum gravity. 

We first draw some general conclusions about the quantity we wish to calculate. 
We know that, by definition, the Hartle-Hawking vacuum respects the symmetries of the 
spacetime and that the metric given by Eq.([3]) is symmetric under the transformations 

T^-T (f)^-(j) e^TT-e, (44) 

and under the translational symmetries 

T ^ T + To + 00. (45) 

We note, however, that unlike the spherically symmetric case, we do not have a 
translational symmetry in the ^-direction, but only the discrete refiective symmetry 
in the equatorial plane. As such, we have the possibility of non-zero off-diagonal 
components T^g and T^r- All other off-diagonal components are identically zero for 
the Hartle-Hawking vacuum. For a massless scalar field propagating in a Ricci-fiat 
background, we have the following expression for the classical stress tensor 

T\ = (1 - 20^7"V;c</';b + (2e - l)6\g'''ip.,ip.,d - 2e^7"V¥';c6 + 2^6\ipnip (46) 

It will prove convenient to write this tensor as 

T\ = [fM^ixMx'))] = lim fMifixMx')) (47) 
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where f^b = T'"'b{x, x') is a differential operator which may be defined in any way provided 
it gives (46) in the coincidence hmit. We shall adopt the following definition, 

f\ = (1 - 20g''''VbV,> + (2^ - '^)6\g'^'V,Vd' - 2^g-'V,Vb + 2e^^V,V^ (48) 

where ^f"*' are the bivectors of parallel transport. Using point separation, the regularized, 
unrenormalized quantum stress tensor in the Hartle-Hawking vacuum state arises when 
this differential operator acts on the thermal Euclidean Green's function (11), 

{f\)unren = (1 - 20^?"^'^;,,, + (2^ - |)5^(7^"'G;ed' - 2^^?"^^;,, + 2^6\DG . (49) 

This is clearly divergent in the coincident limit and must be renormalized prior to 
relaxing the regularization. To achieve this we note that this Green's function is of 
Hadamard form and may be written as 



G{x, x') 



1 



Ai/2 



cr 



+ V\n{a) + W 



(50) 



where a{x,x') is half the square of the geodesic distance between x and x' and 

A = -g-'/\x)det{a,ab')g-'^\x') (51) 



is the Van Vleck-Morette determinant. V{x,x') and W{x,x') are regular symmetric 
biscalar functions with V{x,x) being geometrical (local) while the determination of 
W{x,x') requires (global) boundary conditions. Thus V{x,x) captures the (local) 
ultraviolet divergences of the quantum theory while W{x,x') captures the state 
dependence. 

The divergent part of the stress tensor is then the result of the operator ( 48 ) acting 
on the singular part of the Green's function, i.e. 



div 



-J b 



Ai/2 



cr 



V \n{a] 



(52) 



The renormalized stress tensor then essentially has to be given by the coincident limit 
of the differential operator (48) acting on the regular part 

^ [r\W{x,x')]. 



(53) 



Rather than calculating W{x,x') itself, we calculate (53) by subtracting Eq.(52) from 
(49) and then take thing coincidence. There is one caveat, however, since the application 
of the differential operator f";, to W{x, x') results in a stress tensor which is not conserved 
since W{x,x') is not a solution to the homogeneous wave equation but rather an 
inhomogeneous one where the source is geometrical. It is easy to redeem this situation 
by adding a regular geometrical term to the definition of the stress tensor 

{T\) 



{T\)r 



47r2 



vi6\ 



(54) 



where for a massless scalar field in a Ricci-fiat background, the geometrical scalar vi 
has the very simple form 
1 



(55) 
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On the horizon of our spacetime this reduces to vi{2M) = /€^/(60M^). 

In this paper, we are concerned with the calculation of the stress tensor on the 
horizon only. In this case, it is most convenient to split in the radial direction, taking 
the inner point (say x'), to lie on the horizon. Where derivatives of the Green's function 
are concerned, we must first perform the derivatives before we take partial coincidence 
limits. The bivectors of parallel transport have been calculated for radial separation in 
reference [T] for the Schwarzschild spacetime and are trivially extended to the cosmic 
string case by multiplying the azimuthal component by an factor. We list them here 
for completeness: 

g,,, = (1 - 2M/r)i/2(^ _ 2M/r')'/' grr' = (1 - 2M/r)-^/^{l - 2M/r')"'/' 
gee' = rr' g^^i = rr' sin^ Q. 



6.1. The Derivatives of the Green's Function 



Examination of (49) reveals that there is essentially two types of terms we need to 
evaluate: those of the form g"''^' G-c'b and those of the form g°''^G-cb- 

Considering the former case first, such terms involve a covariant derivative at each 
spacetime point but since G(x, x') is a scalar at both x and x', we are in fact only dealing 
with partial derivatives. For the angular terms g^^'G-eg/ and g'^'^' G .^^^i ^ it is clear from 
the asymptotic forms ( 14 ) that taking x' to lie on the horizon means that only the n = 
terms will contribute in the limit 77' — )■ 1. Therefore, we can differentiate directly Eq.(34) 
which is straight forward. Taking partial coincidence limits and expanding about the 
horizon, we obtain 



and 

[g'^'^'G.,^ 



+ 



^7r2M2(r/ - 1)2 167r2M2(r/ - 1) 327r2M2 
K2(a2 _ 1) . (11^2 + 1) + 30^2 cos2 e 



14407r2M2 



a^sin^^ 



+ 0(r^-l) 



(56) 



87r2M2(r^- 1)2 167r2M2(r^ - 1) 
(a2-l)(a2 + l) 



327r2M2 
+ 0(r7-l). 



+ 480vr2M2^'" a^'^in^; + 0(r^ - 1), (57) 
where we have adopted square bracket notation [..]^ to indicate that we have taken the 
partial coincidence limit (r — )■ r', rj' ^ 1,6 ^ 6' ^ (f)'). 

For g'^'^' G-rr' and g^"^' G-rr'i we must differentiate the full Green's function given by 
Eqs.(ll) and (15) before we can take x' to lie on the horizon. Considering g'^^'G-rr' first, 

00 00 
J2 n^e^"'^^---') Yl 



we have 

9^^ G-tt' 



1 



1 



l=\m\ 



2M/r)V2(i_2M/r')i/2 §7^2 

r(A 
r(A 



fr,\ I i\r(A+|m|/a + l) /)ND-I™I/"^ a'\ i '\ /c-on 

\ I \ — \m\/a + 1) 
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Trivially, the n = term will vanish. Moreover, using the asymptotic forms Eq.(14), it 
is clear that taking x' to the horizon only the terms corresponding n = ±1 survive. We 
now have, in the partial coincidence limit, 

inrr'r 1 (^ + 1)^/^ V-roA I .x r(A + H/« + i) 

iy "^-^rr'lr iQ^2 2^^ (t] - iV/^ ^ ^ ^ ^ ^ T ( A - Iml /« + 1) 

(59) 

We will denote the double sum by 

F(,,cos^)^ f:(2A + l)^^^^^^^[Pr''"''"(cos^)]V(r^), (60) 

m=-oo/=|m| III/ ) 

and return to its evaluation later; this will involve both numerical and analytic parts. 
We follow a similar argument to obtain an expression for g'^'^'G-rr'- We have 

oo oo 

g'''G.rr' = {l-2M/ry/\l-2M/r'Y/^£^ ^ gm«{r-r') J2 e'^^4>-4>') 

71=— oo m=— oo 

Splitting up the = and n 7^ terms and writing in terms of radial variable r/, we 
have 

^rr'r = iv - iy/\v' - ly/' ^ ^ ... r(A+|m|/a + l) 

^ (r7 + l)i/2(Y + 1)1/2 87r2aM3 ^ ^ ^ ^ r(A - |m|/a + 1) 

^ ' ' ^ ' ' m=-oo l=\m\ \ I 1/ / 



xp-l"^l/"rcos^)p-l'"l/"(cos^') f dPA(r70 dgA(r/) ^^n^.-.') dpnA(V) dgnA(r?) 

^ \ dri' dr] dr]' dr] 

' (62) 

Now using the fact that 

(^/2_^)i/2di|(r/) ^pi^^^^^Q as r7'->l (63) 

imphes that the n = term vanishes. Also, using the asymptotic forms Eq.(14), we 
have for n 7^ 0, 

dr]' 2 

which will vanish for all but the n = ±1 terms in the 77' — )■ 1 limit. Therefore, the only 
contribution to g^^'Orr' comes from the n = ±1 terms, and in the partial coincidence 
limit we have 

'^^'J'' 47r2M2a2i/2(r] + 1)1/2 

oo oo -J- 

XA- |m|/a + 1)' ^ dr) 



ATrm^a 2i/2(ry + 1)1/2 
X E E (2A + 1) + 1"^'/" + [p-H/-(eos g)]^%M 

m=-ooi=|m| 

(// — 1)^/^ dP(r7,cosi 
47r2M2a2V2(ry + 1)1/2 



(65) 
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Finally, we require g^"^ G-ye- Similar calculations to those above show that only the 
n = ±1 term survives when we differentiate with respect to r' and then take this point 
to lie on the horizon, yielding 

9^^ G-r'0 



47r2Ma2i/2(^ + 1)1/2 

d V- V-,., ,,r(A+|m|/« + l)^ 

!L Z^(2A + l)— — — — ^P,' (cos0)P^' {cos9)qUri)- (66) 

m=-col=\m\ VII/ / 

In order to take the 6^6' limit, we make use of the fact that 



-oo 



l=\m\ 



(=|m| 

Therefore, in the partial coincidence limit, we have 



K 



2 



{ri - 1)1/2 Q 



= 8.r2Ma2V2(, + 1)1/2 (68) 

We now consider derivatives of the form g"''^G-cb- For such terms, we have two 
derivatives at the same spacetime point, which will involve the Christoffel symbols, 

G-ab = G^ab — F'^bG ,c, (69) 

where the Christoffel symbols are the same as those of Schwarzschild spacetime except 

r^_^ = -a2(r -2M)sin2^ TJ^ = -a^ sin cos ^ . (70) 

In spite of this minor added complication, things are significantly easier in this case since 
only the n = term will contribute for each derivative. This is because we can always 
choose the differentiation to act on the outer point, this will not affect the asymptotic 
structure at the other spacetime point and so taking 77' — )• 1 and using the asymptotic 
forms (14), it is clear that only the n = term will contribute in the limit. Now, the 
derivatives in (69) can be obtained by direct differentiation of (34). Performing the 
derivatives and using the appropriate Christoffel symbols, followed by taking the partial 
coincidence limits and expanding about the horizon, we obtain 



87r^M^{r] - 1)2 87T^M^{ri - 1) 327r2M2 



14407r2M2 a4 gij^4 q 



+ 0(r^-l), (71) 



^ Snm^iri - 1)2 87r2M2(r/ - 1) 327r2M2 



,^(^l^i)(i]f^,0(,-l), (72, 

144U7r^ikZ ^ a* sm 6 
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k2 (a2-l)(49a2-l) 



14407r2M2 a4 51^4 q 



+ 0(r/-l), (73) 



/t^ ^ ft:2 ^2 (a2_i)(«2 + i) 



87r2M2(r7 - 1)2 327r2M2 ASOtt^M^ a^sin^O 
k"^ (a^-l)cos^6' ^, 
Aqu'^M'^ sm y 

, /t2 (a2 - 1) (11^2 + 1) cos g /t2 (a2_i)cosg _ 

7207r2M3 a^sin^^ 967r2M3 a2sm3^ ^ ^' 

(75) 

^j.^. Candelas Method for Obtaining a Series Solution for F{ri, cos 9) 

In order to proceed with the evaluation of the stress tensor on the horizon, we require a 



series solution for the function we have called F{r], cos 9) (60) about the horizon 77 = 1. 
Fortunately, Candelas [T] has obtained the series solution in the Schwarzschild case and 
the method described here follows an almost identical path, so we refer the reader to the 
Candelas paper for more details. Following Candelas, we may capture the asymptotic 
behaviour of qix{ri) as — ?■ 1 by writing 

qixiv) = -V2Q\{r]) + PxPixiv) - f{v)Qxiv) + Qixiv) (76) 

where the (3\ are constants depending on / and m which can be evaluated numerically 
and the function f{ri) and the f3x are chosen such that giA does not contribute as — )■ 1. 
The /3\ coefficients are global in nature and depend on the boundary conditions imposed, 
as we shall see in Appendix C. Further following Candelas, we may take 

V2 r ((e + 1)^-16) ,, 

fiv) = ^j^ (e- 1)3/2 (77) 
so, in particular, have 

fiv) = iv- 1)'/' + UV- + Oiv - ly^' . (78) 

A Frobenius analysis on this differential equation about 77 = 1 now shows that qi\{ri) 
satisfies the asymptotic condition 

~ ihiv - ly^'Hv - 1) as r/ ^ 1, (79) 

which is uniformly valid in / and m. As a consequence, Candelas [1] has shown that for 
this choice of /(r/) and the coefficients (3x, the remainder term qi\{rj) does not contribute 
in the limit — )■ 1; although Candelas's analysis was for A integer, the argument applies 
immediately to non-integer A. The contribution coming from qi\{r]) + PiPixif]) can now 



be calculated to leading order in {rj — 1), using the asymptotic forms (14) for pix and 
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the fact that qw does not contribute in the hmit, as 

oo / 



(2A + i)r^_|^|/^ + ia^A 



-\m\/a 



COS I 



i2r~ 



i=0 m=-i 



where 



(A - \m\/a + 1) 
[r] - iy/^A{cose) + 0{7] - 1) 



[qixiv) + l^xPixiv)] 



-1^3/2 



(80) 



A cos^ = V V 2A + 1 ^ , \, 
^ ' ^ ^r A- m/a + ir^ 

This sum is highly convergent and easily evaluated numerically. As we show in more 
detail in Appendix C, the (i\ coefficients are determined by equating the asymptotic 
expansions for two equivalent representations of the function q\\{j]\ one being an 
integral representation obtained from the Wronskian condition and the other being the 
form (76). Returning now to the evaluation of F{ri,cos6), 

oo I 

F{ri, cos 9] 



COS I 



'/3a. (81) 



1=0 m=-l ^ I 1/ < / 

1=0 m=-l ^ I 1/ "I" J 



-V2{r] -1)/ + qix{ri) 



(82) 



For the second and third terms in the expression above, we can use the generalized 
Heine identity (32) and its derivative respectively. As a series expansion in [rj — 1), we 
obtain 



F(?7, cos( 



2a 



{r, - 1)3/2 2(77 - 1) 
(1 + 10^2 - iia^ 



^-^ + („-l)"=(^(cos 



1) 



a [a 

V - — 

16 6a sin"* 9 



+ 0{ri 



\3/2 



(83) 



90a3 sin4 q 

We can now calculate the asymptotic series expansions for g'^'^'G-rr', g"'G-rr' and 



g'^'^'G-r'0. Employing Eq.(83) in expressions (59), (65) and (68) yields 



K 



{a' - 1) 



Sn'^M^ir] - 1)2 647r2M2 967r2M2 ^2 sin^ 9 



[a 



1)(11«=^ + 1) 



14407r2M2 a^sm^9 
3k^ 



+ 



K 



+ 



167r2M2a 



A(cos0) + 0(r/- 1), (84) 
^2 (a2 _ 1) 



+ 



87r2M2(r/- 1)2 87r2M2(r/ - 1) 647r2M2 967r2M2 ^2 sin^ 5 



(a2- l)(lla2 + i) 



14407r2M2 
W^'G.Mr = 0{r,-l). 



a^sin^^ 



+ 



167r2M2a 



A(cos^) + 0(7]- 1), (85) 

(86) 
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6.3. Renormalized Stress Tensor 



We have now calculated all the derivatives of the Green's function we require in order 
to calculate the unrenormalized components of the stress tensor. The term g'^^'G-ab' 
appearing in the stress tensor (49) is found by summing the expressions (84), (85), (56) 
and (57), yielding 



ab'W 



y4(cOS( 



K 



247r2M2 



sin*^ 



O(r^-l). (87) 



Moreover, the term DG appearing in (49) vanishes since, by definition, G{x,x') for 
a massless, Ricci-flat scalar field is a solution to the inhomogeneous wave-equation 
□G = —S. There is, however, a contribution to the divergent subtraction terms (52) 
coming from nGdiv since this is non-zero. Combining these results, we arrive at the 
following expression for the unrenormalized stress tensor components near the horizon: 



K 



,^2 



2e 



1 



87r2M2 {r] - 1) 



+ (14^ - i; 



K 



647r2M2 



[a 



1){1W + 1] 



14407r2M2 



3k^ 



sin^ 9 

(a^ - 11 



2'247r2M2 a2sin^^ 



K 



1 



K 



I) 

6n67r2M2 a2 sin^^' 

.2 



^7r2M2 {r] - 1)2 

,^2 



?7r2M2 {ri 



+ (2^-3); 



K 



;7r2M2a 



y4(cOS( 



K 



1) ' ' ^ "'647r2M2 
(a2-i)(iia2 + l) 



' 14407r2M2 



a4 sin^ Q 



K 



1) 



2'247r2M2 d^mi^Q 
1 



87r2M2 {^ - 1)2 



+ (2e 



87r2M2a 
{o? - 11 



- (1 - 2^) 



-A(cos I 



1^ 

6^167r2M2 a2 sin^^' 

1 



167r2M2 {^ - 1) 



(1-40: 



327r2M2 



14407r2M2 



247r2M2 cfi sin^ Q 

1 



87r2M2 [ri 



1-20 



1 



167r2M2 {ri - T 



sin^ ^ 



+ (1 - 40 



327r2M2 



2'87r2M2a 



A(cos6') + 



(a2_i)(a2 + l) 



4807r2M2 a4 51^4 q 



\a 



1) COs2 d 



6'87r2M2 a2 ^^^^ 
[a^ - + l)cose 



3607r2M 



sin^ ^ 



K 



(a 



1) cosO 



sin'^6' 



(88) 



(89) 



(90) 



(91) 



(92) 



Note that the off-diagonal component is regular and vanishing for minimal coupling. 
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The geometrical subtraction terms are found by obtaining a series expansion for 
the singular part of the Green's function, then applying the differential operator (48) 
and taking partial coincidence limits; this is completely equivalent to, for example, 
Christensen's method |9]. The results are 



{T r) div 

div ■ 
{T^e)div ■ 
div 



K 



Stt^M^ {t] - 1)2 
3k2 1 



+ 



,7r2M2 {t] - 1) 

1 



TT^M^ {t] - 1)2 87r2M2 {t] - 1) 



K 



K 



1 



div 



87r2M2 [r] - 1)2 

/t2 1 

' 87r2M2 [t] - 1)2 
0. 



(1 - 2^) 



167r2M2 (ry - 1) 
1 

167r2M2 ij] - 1) 



/t2 

4'l67r2M2' 



K 



T_ 

20^167r2M2^ 

2 



K 



20' 



167r2M2' 



(93) 

(94) 

(95) 

(96) 
(97) 



Finally, from Eqs.(54), (88)-(92) and (93)-(97), we arrive at the renormalized stress 
tensor for a massless scalar field in the Hartle-Hawking vacuum on the horizon of the 
Schwarzschild black hole threaded by an infinite cosmic string: 



15 ■^^327r2M2 ' ■^87r2M2a 



-(2^ 



K 



K 



2 



A(cos I 



K 



(a2_i)(ii«2^1) 



14407r2M2 



sin^ e 



[a 



2'247r2M2 a2sin^^ 



6^167r2M2a2sin2 0' 



(98) 



15 



'13 
^60 



-0 

■(2e 



-0 



327r2M2 ^87r2M2a 



-A(cos^) 



(a2-l)(lla2 + l) 



14407r2M2 



sin^ e 



[a 



2'247r2M2 a2sin4^ 



1) 



167r2M2 a2 sin^^' 



(99) 



K 



167r2M2 



+ (2^ - I 



K 



A{cos9) 



(a2_ i)(iia2 + l) 



14407r2M2 



K 



(a2 - 1) 



247r2M2 a^sin^e' 

2 



^60 



K 



167r2M2 



+ (2e 



1' 



K 



[a 



1) cos^ 9 



a 



-A{cos9) + 



sin^ 9 

(100) 

(«2_i)(a2^1) 



4807r2M2 sin^ ^ 



K 



6'87r2M2 a2 sin^^' 
2 (a2-l)(lla2 + l)cose 



fa^ - 1) cos^ 



(101) 



(102) 



'3607r2M a^six^'9 ^ iSn^M sin^ 9' 

We note from these expressions that {T'^r)ren = {T^r)ren, which we would have 
expected from the symmetries of the black hole horizon. We also see that {T^e)ren 7^ 
{T'^^)ren but that wc retrieve equality between these components in the a — 1 hmit. 
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0.6r 




(87r^MV)<r"«>,e„ 
(8;r2M-/r')<r%>„„ 



Figure 2. (87r2Af2/K2)(T°6)J?°;*^°" 
as a function of cos 6 ior a ~ 0.95. 

(8;:^MV)<r°i,>re 



for a conformally coupled scalar field, ^ — 1/6, 
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-0.04 





(87r'M'/K')<rr>„„ 
(8;r'MV)<-J''»>„ 
(87r2MV)<r"o>,e„ 
(8;r^*fV)<3"%>re„ 



Figure 3. {Sn'^ IvP / K'^){T''b)'^°"^°" (zoomed in plot) for a conformally coupled scalar 
field, ^ — 1/6, as a function of cos0 for a = 0.95. 



Furthermore, in the conformally coupled case, ^ = 1/6, we reproduce the Candelas [T] 



result in the a — )■ 1 limit. Eqs.(98)-(101 ) also, of course, yield new results for the 



renormalized stress tensor on the horizon of a Schwarzschild black hole for general ^, 
these are trivially obtained from by taking the limit a — )■ 1. 

In Fig. |2] and Fig. [6] we have plotted the components of the renormalized stress 
tensor for the massless, conformally coupled scalar field for a = 0.95 and a = 0.6, 
respectively. We have further plotted the minimally coupled case in Fig. |4] and Fig. [7] 
for a = 0.95 and a = 0.6, respectively. We see that for small azimuthal deficits, the 
values only move away from the constant Schwarzschild values very close to the poles 
where the curvature singularity dominates. For this reason, we have included magnified 
plots Fig. |3]and Fig. [sjfor the a = 0.95 case. For larger deficits (a = 0.6 for example), 
we observe that the effects of the string are apparent even relatively close to the equator. 
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(8rMV)<7'°i,> 
1.0 




(87r2MV)<r^>re 

(87r2MV)<r"o>,e 
(8;i2MV)<7'%>r, 



Figure 4. {8Tr'^AP/K'^){T"-f,)ren^°" a minimally coupled scalar field, ^ = 0, as a 
function of cos^ for a ~ 0.95. 
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Figure 5. {Sn'^M'^ / k^){T<',,)^°"''°'^ 



(zoomed in plot) for a minimally coupled scalar 



field, ^ = 0, as a function of cos 6* for a = 0.95. 
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Figure 6. (87r2MVK2)(T%)J?°;*^°" 
as a function of cos 6* for a — 0.6. 



for a conformally coupled scalar field, ^ = 1/6, 
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Figure 7. {Stt'^ AP / k'^) {T" t)';:";;^'"'"' for a minimally coupled scalar field, ^ = 0, as a 
function of cos^ for a ~ 0.6. 



7. The Conservation Equations: The Minimally Coupled Case 

In this section, we give an example of how the stress tensor we have calculated may be 
shown to satisfy the conservation equations. For simplicity, we consider the minimally 
coupled case ^ = 0. The conservation equations are 

Va(T'^fe)ren = da{T"'b)ren + ^lA^'^ b) ren — ^'ab(^°' c) ren = 0. (103) 

In particular taking b = 6 we have 

^a{T"'e)ren = de^T^ e)ren + dr{T^ e)ren H {T^9)ren + {r — 2M){T^r)ren + 

r 

cot e{f\)ren ' COt e{f^^)ren = 0. (104) 



Now 



(r - 2M) {f\)ren = " {f'e)ren ^0 as r ^ 2M (105) 
r 



for the minimally coupled case. However, the derivative dr{T^e)ren is non-zero on the 
horizon, as we shall see below. For = 0, the classical expression for this term is given 
by 

dr{f'~e)ren = (9^ (/''v2,r)V5,e + /""V.rV^.re (106) 

and in the regularized version terms like (p^r'P,9 are reexpressed with a radial derivative 
of the Green's function at one spacetime point and a 6 derivative at the other spacetime 



point. We consider the first term on the right-hand side of Eq.(106), taking the radial 
derivatives with respect to the outer point and taking the inner point to lie on the 
horizon, then only the n = term survives and we can differentiate our Heine identity 



(|32|) directly 

2Ti'^M'^a drj \{ri + 1) drjdO' \sin^sin^'sinhx(cosh(x/Q;) — 1) 



^/rr N d f {r]-l) f sinh(x/a) \\ (^n7\ 
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where 



coshx 



77 — cos 6 cos 6' 



sin 6 sin 6' 
and we have taken t ^ r' and — 
taking 6 ^ 6' , we have 

[(9r(/>,r)</',e]r ^ 



(108) 

Finally, performing the differentiation and 



[a 



l)(lla2 ^ i)cos^ 



3607r2M2 



sin^ ^ 



+ 0(r7- 1)^/^(109) 



In taking the 6^6' limit, we have made use of Eq.(67). 



For the second term on the right-hand side of Eq.( |106[ ), we have 

d d d 



00 



00 00 

n=~oo m=—oo l=\m\ 



X 



r(A+ |m|/a + l) 



p-l'"l/"(cos^)p-l™l/"(cos^')XnA(r,r'; 



;iio) 



r(A- |m|/a + l)' 

where we have taken the partial limits — )■ 0', r — t- t'. Taking the 'primed' coordinates 
to correspond the inner spacetime point, we can see from the asymptotic forms (14) that 
taking r' to lie on the horizon eliminates all but the n = ±1 modes. The contribution 
coming from these modes is 

^2 (^2 _ 1)1/2 g Q 



K 



8\/2 7r2M2a (r/ + 1) dri 09 



F{ri, cos 6) 



111 



where F(?7,cos6') is given by Eq.(82) and in going to the limit 6 ^ 6' we have again 
made use of Eq.(67). Now, using our series solution (83) for F(?7,cos^), we arrive at 



K 



d 



K 



[a 



A{cos9) + 
- 1) cos 9 



(a 



+ l)cos^ 



487r2M2 



sin'^6' 



3607r2M2 



sin^ 6 



Combining these results, we have 



[dr{T'^ e)ren\ 



K 



d 



167r2M2a 89 



A{cos9) 



[a 



2 - 1) cos^ 



487r2M2 



sin"^^ 



(112) 



;ii3) 



There are no subtraction terms for this quantity so that this expression represents the 
renormalized expectation value dr{Tg)ren- Substituting Eqs.(113), (100) and (101), for 
,^ = 0, into Eq.(104), it is readily seen that 

= 



as required. One can similarly, show that Va{T"'r)r 
are trivially satisfied. 



(114) 

0. The equations for r and 



8. Conclusions 



We have considered a massless, arbitrarily coupled quantum scalar field on a 
Schwarzschild black hole threaded by an infinite straight cosmic string. We have 
obtained simple expressions for the both the renormalized vacuum polarization and 
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the renormalized stress tensor on the black hole horizon for this field in the Hartle- 
Hawking vacuum state, which is a two-fold generalization of the Candelas [1] result to 
encompass an arbitrary coupled field and a black hole threaded by a cosmic string. In 
order to do this, we have derived a very useful summation formula involving non-integer 
Legendre functions, which we have called a generalized Heine Identity given its analogy 
with the standard Heine Identity involving Legendre functions of integer degree. The 
derivation of this formula involved equating equivalent expressions for the same Green's 
function on a conveniently chosen space-time. Furthermore, the derivation did not rely 
on specific properties of the Legendre functions. This idea can be used to derive many 
useful summation formulae and addition theorems for non-integer Legendre functions, 
some of which are presented in Appendix B. Elsewhere [11], we have used this idea to 
obtain an expression for the vacuum polarization on the entire exterior region of the 
space-time. We believe this framework will prove useful in the other axially symmetric 
black hole calculations, most importantly, the Kerr-Newman black hole. 



Appendix A: Normalization Condition 



In this Appendix, we shall prove the normalization condition Eq.(lO) for the non- 
integer Legendre functions that arise on separation of the wave-equation on cosmic 
string spacetimes. We begin by noting the following relationship between the Legendre 
function and the Gegenbauer polynomial [T2j, 

2V2-'^r(2i^)r(n -fl) 



p'/'"" {x^ 



(A.i; 



where C^{z) is the Gegenbauer polynomial and n is an integer. We also have a 
normalization for these polynomials [12] given by 



2i-2^'vrr(n + 2i/) 
(n + u)r(n + l)r(u) 



(A.2) 



Choosing n = I — \m\ and u = l/2 + |m|/ain Eq.( |A.l ) gives us the appropriate 
Legendre function in terms of the Gegenbauer polynomials, and with this choice of n 
and u, Eq.(A.2) becomes 



p-\m\/a , 



cos^)P^7i/:i^l/Jcos^)d(cose) 



2-4 


'^l/"7rr(l + 2|m| 


/af 2 r(A - 


m\/a + 1) 


r(i + 


m|/a)2r(l/2 + 


m|/a)2(2A + l)r(A + 


m\/a + 1) 



5iv 



We now use the so-called doubling formula for the Gamma function [16 

22|m|/, 



r 1 + 2 m /a 



TT 



-r(l/2+ |m|/a)r(l + \m\/a) 



which implies that 

2-4|m|/a^PQ ^2|m|/a)2 



r(l + |m|/a)2r(l/2 + |m|/«) 



1. 



(A.3) 



(A.4) 



(A.5) 
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Finally, substituting this into Eq.(A.3) gives the desired result, viz. 



2 r(A-|m|/a + l 



(2A + l)r(A+|m|/a + r 
Appendix B: Summation Formulae 



6w. (A.6) 



In this Appendix, we shall derive some interesting summation formulae for non-integer 
Legendre functions based on considering equivalent forms of the 3D Green's function on 
flat space threaded by a cosmic string 

ds^ = dr^ + rM^^ + aV^ sin^ Odcf^ . (B.7) 
The 3D Green's function for Laplace's equation satisfying 

V2G3(x, x') = -(73(x)-^/^5(x - x') (B.8) 
can be written in spherical polar coordinates as the mode sum: 

r V p^™(*-<^') V r(A + |m|/a + i) 

' '^>- Aira ^ ^ r(A - \m\/a + 1) 

m=-oo l=\m\ ^11/ ' 

P;l"^l/"(cos^)P;l'"l/"(cos0')(^<)\^>)"^^^'^- (B.9) 

Note that this corresponds to the ra = term of the 4D Green function in the 
corresponding ultrastatic {gu = — 1) space-time, but that the corresponding statement is 
not true in the merely static Schwarzschild cosmic string space-time. We can also write 
this same Green's function in several other forms and in several coordinate systems, for 
example, performing the mode decomposition in cylindrical polar coordinates, p = r sin 6 
and z = r cos 6, we have 

G3(x,xO = ^ J2 e™(<^-^') / coskiz-z')Ii„,yaikp^)Kim\/aikpy)dk. (B.IO) 

m=— oo ^ 



The integral over k may further be performed to yield 

m=— OO ^ ' 



Another form of this Green's function is obtained from General Axi-Symmetric Potential 
Theory [HI [18]: 



m=— oo 



/o [(2-2')^ + P^ + P'^-2pp'cos^](H/«+i/2)^^- (^-^^^ 

These expressions are valid for all a G M+; equating any one of these expressions with 
another gives us a formula for the functions involved. Some of these results are well 
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known and can be found in volumes such as Gradsteyn & Rhyzik fTB], nevertheless, 
we have non-obvious, important results, for example, equating (B.9) with (B.ll) and 
taking r — )■ r', we obtain 



E 

l=\m\ 



\m\/a 



r(A- |m|/a + 1) 



cos^O 



TT sm 



in^sin^')V2^l"l/"''/' V sin ^ sin ^' 



1 — cos cos 0' 



(B.13) 



Also, equating Eq.(B.9) with Eq.(B.12) we arrive at the following formula 



l=\m\ 



:d^. 



[{Z - Z'f + p2 + p/2 _ 2pp' cos ^](l™l/«+l/2) ^ ^ ' ^^-^^^ 

SO that the integral expression on the RHS is a type of generating function for these 
non-integer Legendre functions. 

There are also a set of formulae that hold for a more restrictive set of a values. 
Linet [19] has shown that, when a > 1/2, the Green's function may be written as 



G3(X,X') 



1 1 

4^ ^ 8 



(B.15) 



where 



a = [p2 + p'2 _^ _ ^/^2 _ 2pp' cos - 
= [p2 + p'2 + _ ^/)2 ^ 2pp' cosh m] 



sin(\E' — Ti/a) 



sin(\E' + tt/q;) 



and 

' cosh(M/a) — cos(\E' — tt/q;) cosh(u/Q;) — cos(\E' + tt/q;) 

Equating this to ( |B.9[ ) for a > 1/2 and taking r — )■ r', we find 



(B.16) 
(B.17) 



m=— oo Z=|m| 



r(A + \m\/a + l) 
f(A 



P. 



-|m|/o 



+ 



\m\/a + 1 
1 1 



A 



COS^)P. 



— |m|/a 



A 



COS 0') 



where 



[2(l-cosa7)]V2 27ra (2sin^sin^')^/^ Jo [cosh ^ + cosh 

cos a7 = cos 9 cos 9' + sin 9 sin 6'' cos a{(f> — 0') 

_ ^ 1 — cos 9 cos 6*' 
cosh 4 = 



(B.18) 



sm sm 9 

This formula has the benefit of relating a 3D mode-sum expression to an expression that 
isolates the 3D Hadamard singularity structure (the first term on the RHS) together 
with a regular term that depends on the boundary conditions. We have used this result 
to prove the regularity of the mode-sum that arises in the calculation of the vacuum 
polarization on the exterior region of the spacetime [TT]. 



Schwarzs child Black Hole Threaded by a Cosmic String 25 

Another form of the Green's function can be given in terms of the half-integer 

Legendre functions by considering the Green's function in toroidal coordinates, which 

are related to the cartesian coordinates by 

sinh u cos sinh u sin sin n 

X = — y = — z = — , B.20 

cosh /i — cos rj cosh n — cos rj cosh n — cos rj 

where the ranges of the coordinates are 0</i<oo,0<?7<27r,0<(/)<27r. The 
mode-sum form of the Green's function obtained by separating in these coordinates is 

Gs(^,^') = —- y e*-(<^-*') V e'"(''-'''H(cosh/i- cos r/) (cosh /i'- cos r/')]'^' 



P„-^X(cosh/i<)g;i7y,"(cosh/i>). (B.21) 



T{n-\m\/a + 1/2) ""^/^ ^ '-<^^n-i/ 



Equating this to (B.ll) gives us the following addition theorem 

^ ,.«,,-,',,.iH./«Ei!i±M^±M 



n=— oo 



where 



r(n- |m|/a + l/2) ""^/^ 
(sinh ;z sinh /.OV^^I"^!/"-^/^^^^- ^^'^^^ 

^ 2pp' ^ ^ 

which are related to the toroidal coordinates by 

p= sinr^ _ 

cosh /i — cos rj cosh /i — cos rj 

The final coordinate system we shall investigate is prolate spheroidal coordinates 
which may trivially extended to the oblate spheroidal case. The relationship between 
prolate spheroidal coordinates and cartesian coordinates is given by 

a; = sinh (T sin 6* cos ?/ = sinh a sin ^ sin z = cosh cr cos ^. (B.25) 

where the ranges of the coordinates are O<cr<oo,O<^^<7r,O<0<27r. The mode 
form of the Green's function obtained by separating in these coordinates is 

^3lx,xj-^^^ e 2^[ZA + L)e r(A-|m|/« + l)2 

m=-OD l=\m\ \ I II I 

p-l'"l/"(cos0)P^"l'"l/"(cos0')^A~'™'^"(cosha<)g^'"l/°(cosha>) (B.26) 

We shall write down only one possible summation formula here, by equating the Green's 
function expression above to (B.ll), we obtain a summation formula for a product of 
four Legendre functions: 

P,-l'"l/"(cosha<)Q-l'"l/"(cosh(7>) = - ^ ^g^?) 

TT V sinh cr sinh & sin Q sin Q' 
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where 

cosh^ a + cosh^ a' — s\v? 9 — sin^ 9' — 2 cosh a cosh a' cos 9 cos 9' ^g-j 
2 sinh cr sinh a' sin 6' sin 6'' 
There are many more summation and addition formulae that can be derived in this way. 
Those presented here were either of particular use to us elsewhere or are to the best of 
our knowledge previously unpublished. 



Appendix C: Determining the j3x Coefficients 



We briefly describe the procedure for determining the (3\ coefficients occurring in the 
definition of A{cos9). The method described is completely analogous to the Candelas 
|1] method except that we replace / by A = / — \m\ + |m|/a everywhere it occurs. We 
begin with the Wronskian relation between the n = ±1 radial functions of the first and 
second kind, 

(C.29) 



Dividing across by {pixY and integrating we obtain the following integral expression for 



Qixiv) = 2pia(^) 



(C.30) 



(e-i)biA(O]^ 

where we have used the fact that qixijf) — )■ as ?7 — )• oo. One can easily verify that the 
series solution to about ^ = 1 is 

pixii) = {i- If" + \K\ + m - ^f" + o(e - 1)'/' (C.31) 

and therefore we have 

1 1 



+ o(e-i: 



(C.32) 



We now add and subtract the first two terms on the right-hand side above in the integral 
expression (C.30), 



qix{ri) = 2pix{r]) 

+ 2piA(r?) 



1 



1 



(e-l)V[PlA(eP (e-1) 

de f I A(A + 1 



+ 



A(A + i: 



- 1) - 1 



(C.33) 



It is clear from Eqs.(C.31) and (C.32) that the first integral is convergent as r/ — )■ 1, and 
the second integral can be evaluated explicitly. We thus obtain the following asymptotic 
expansion of ^^(r]) as ?7 — )■ 1, 

^ +1[1 + A(A + 1)](7^- 1)1/2 ln(7^-l) 



Qixiv) 



(r/- 1)1/2 2 1 

{2/, + |A(A + 1) - i[l + A(A + 1)] ln2}(r^ - if/' + 0((r/ - 1)'/Hn{r^ - 1)) , (C.34) 
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where 

/ 1 1 A(A + 1)' 



I 



A 



Recall that we also had an alternate expression for qix{r]) (76) which we expand 
about rj = 1 to get 

^i^(^) = _\y/2 + + + 1)K^ - 1)'/' - 1) 

+ {[l + A(A + l)][vl/(A + l)+7-^ln2-l] + i + /3,}(r/-l)V2 

where "^{z) = r'[z)/r{z) and 7 is Euler's constant. 



Table CI. The coefficients /3a for a = 0.95. 





m= 


1 


2 


3 


4 


5 




0.119377314 












/=1 


0.002779252 


0.002318785 










1=2 


0.000144497 


0.000127179 


0.000112214 








1=3 


0.000018348 


0.000016749 


0.000015310 


0.000014014 






1=4 


0.000003975 


0.000003705 


0.000003456 


0.000003226 


0.000003013 




1=5 


0.000001181 


0.000001115 


0.000001053 


0.000000995 


0.000000941 


0.000000890 



Table C2. The coefficients for a = 0.6. 





m= 


1 


2 


3 


4 


5 


/=0 


0.119377314 












/=1 


0.002779252 


0.000344996 










1=2 


0.000144497 


0.000033858 


0.000010533 








1=3 


0.000018348 


0.000006342 


0.000002579 


0.000001181 






1=4 


0.000003975 


0.000001723 


0.000000828 


0.000000431 


0.000000239 




1=5 


0.000001181 


0.000000592 


0.000000319 


0.000000182 


0.000000109 


0.000000068 



Finally, comparing our two equivalent expressions for qixirf) given by (C.34) and 



(C.36) yields an expression for fix that is numerically tractable: 
fix = 2Ix - [1 + A(A + 1)] [^(A + 1) + 7 - I 



1 

2" 



(C.37) 



That the derivation of the /3x coefficients depends on the integral (C.30) underlines 
the fact these coefficients are truly global in nature since the dependence on the boundary 
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conditions is explicit here. If we considered the case where we place the black hole in a 
spherical box, then the upper limit of the integration above would be the radius of the 
box and we would obtain different f3\ coefficients. 

The integral Ix is most effectively calculated by splitting the range of integration 
into the regions 1 < ^ < 2 and 2 < ^ < oo. In the former range, one can use a series 
solution to pix{ri) in order to cancel the divergences in the integrand near 77 = 1, and 
the integral can now be performed without difficulty. In the latter range, the first term 
cuts off exponentially as 77 — )■ 00 since pixirf) — )• 00 as 77 — )■ 00. The remaining part of 
the integral can be performed explicitly. We have tabulated above some Px coefficients 
for a = 0.95 and a = 0.6 for / up to / = 5. We note that for m = 0, we have A = / an 
integer, and we retrieve the Candelas [l\ /3i values. 



Acknowledgements 

PT is supported by the Irish Research Council for Science, Engineering and Technology, 
funded by the National Development Plan. 

The authors would like to thank Cormac Breen for many helpful conversations. 



References 



P. Candelas, Phys. Rev. D, 21:2185, 1980. 

P. Candelas and K.W. Howard, Phys. Rev. D, 29:1618, 1984. 

P.R. Anderson, Phys. Rev. D, 39:3785, 1989. 

B.P. Jensen and A.C. Ottewill, Phys. Rev. D, 39:1130, 1989. 

J.G. McLaughlin, B.P. Jensen and A.C. Ottewill, Phys. Rev. D, 45:3002, 1992. 

J.G. McLaughlin, B.P. Jensen and A.C. Ottewill, Phys. Rev. D, 51:5676, 1995. 

P.R. Anderson, Phys. Rev. D, 41:1152, 1990. 

E. Winstanley and P.M. Young, Phys. Rev. D, 77, 2008. 
S.M. Christensen, Phys. Rev. D, 14:2490, 1976. 

P.C.W. Davies and V. Sahni, Class. Quantum. Grav, 5:1-17, 1988. 

A.C. Ottewill and P. Taylor, arXiv:1007:0051[gr-qc], to appear in Phys. Rev. D. 

F. Oberhettinger, A. Erdelyi, W. Magnus and F.G. Tricomi, Higher Transcendental Functions. 
McGraw-Hill, New York, 1953. 

A.G. Smith, Formation and Evolution of Cosmic Strings, eds G. W. Gibbons, S. W. Hawking 

and T. Vaschaspati, (Cambridge University Press), 1990. 
Y. Decanini and A. Folacci, Phys. Rev. D, 78:044025, 2008. 
M.R. Brown and A.C. Ottewill, Phys. Rev. D, 34:1776, 1986. 

I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series and Products. Academic Press, 2000. 

A. Weinstein, Amer. Math. Soc, 63:342, 1948. 

B. Linet, Phys. Lett, 60A(5):395, 1977. 
B. Linet, Phys. Rev. D, 35:536-539, 1987. 



